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Abstract 



We present a measurement of the Standard Model CP violation parameter 
sin 20i (also known as sin 2/3) based on a 10.5 fb _1 data sample collected at the 
T(4S) resonance with the Belle detector at the KEKB asymmetric e + e~ col- 
lider. One neutral B meson is reconstructed in the J/^Ks, ip(2S)Ks, XciKs, 
rj c Ks, J/tpK^ or J/tpir CP-eigenstate decay channel and the flavor of the 
accompanying B meson is identified from its charged particle decay products. 
From the asymmetry in the distribution of the time interval between the two 
P-meson decay points, we determine sin20i = 0.58to;34(stat)taio(syst). 

PACS numbers:11.30.Er,12.15.Hh,13.25.Hw 
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In the Standard Model (SM), CP violation arises from a complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix [0]. In particular, the SM predicts a CP 
violating asymmetry in the time-dependent rates for B® and 1^° decays to a common CP 
eigenstate, fcp, without theoretical ambiguity due to strong interactions 0: 



where T(Bd (B®) — ► fcp) is the decay rate for a Bd (B®) to fcp at a proper time t after 
production, £y is the CP-eigenvalue of fcp, Am d is the mass difference between the two B® 
mass eigenstates, and <pi is one of the three internal angles of the CKM Unitarity Triangle, 



In this Letter, we report a measurement of sin 20i using B^Bd meson pairs produced at 
the T(4S) resonance, where the two mesons remain in a coherent p-wave state until one of 
them decays. The decay of one of the B mesons to a self-tagging state, f tag , i-e. a final state 
that distinguishes between B% and Bf, at time t tag projects the accompanying meson onto 
the opposite 6-flavor at that time; this meson decays to fcp at time tcp- The CP violation 
manifests itself as an asymmetry A(At), where At is the proper time interval At = tcp—tt ag . 

The data sample corresponds to an integrated luminosity of 10.5 ftT 1 collected with the 
Belle detector gj at the KEKB asymmetric e+e" (3.5 on 8 GeV) collider @. At KEKB, the 
T(4S) is produced with a Lorentz boost of = 0.425 along the electron beam direction (z 
direction). Because the B® and Bf mesons are nearly at rest in the T(4S) center of mass 
system (cms), At can be determined from the z distance between the fcp and f tag decay 
vertices, Az = zcp — z tag , as At ~ Az/P^c. 

The Belle detector consists of a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of 1188 aerogel Cerenkov counters (ACC), 128 time- 
of-flight (TOF) scintillation counters, and an electromagnetic calorimeter containing 8736 
CsI(Tl) crystals (ECL) all located inside a 3.4-m-diameter superconducting solenoid that 
generates a 1.5 T magnetic field. The transverse momentum resolution for charged tracks 
is (<y pt /pt) 2 = (0.0019p t ) 2 + (0.0034) 2 , where p t is in GeV/c, and the impact parameter 
resolutions for p = 1 GeV/c tracks at normal incidence are oy^ ~ a z = 55/im. Specific 
ionization (dE/dx) measurements in the CDC (cr dE/dx — 6.9% for minimum ionizing pions), 
TOF flight-time measurements (<Jtof — 95 ps), and the response of the ACC provide K ± 
identification with an efficiency of ~ 85% and a charged pion fake rate of ~ 10% for all 
momenta up to 3.5 GeV/c. Photons are identified as ECL showers that have a minimum 
energy of 20 MeV and are not matched to a charged track. The photon energy resolution is 
(cte/E) 2 = (0.013) 2 + (0.0007/P) 2 + (0.008/ J E 1 / 4 ) 2 , where E is in GeV. Electron identification 
is based on a combination of CDC dE/dx information, the ACC response, and the position 
relative to the extrapolated track, shape and energy deposit of the associated ECL shower. 
The efficiency is greater than 90% and the hadron fake rate is ~ 0.3% for p > 1 GeV/c. An 
iron flux- return yoke outside the solenoid, comprised of 14 layers of 4.7-cm-thick iron plates 
interleaved with a system of resistive plate counters (KLM), provides muon identification 
with an efficiency greater than 90% and a hadron fake rate less than 2% for p > 1 GeV/c. 
The KLM is used in conjunction with the ECL to detect Kl mesons; the angular resolution 
of the Kl direction measurement ranges between 1.5° and 3°. 



A(t) 



r(B d ^f CP )-T(B° d ^f C p) 



£f sin 20i sin Am^i, 
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We reconstruct B® decays to the following CP eigenstates: J/ipK s , ip(2S)K s , XciKs, 
r] c Ks for = —1 and J/ipn , J/ipKi for = +1. The J/ip and ip(2S) mesons are recon- 
structed via their decays to (£ = /i, e). The ij){2S) is also reconstructed via its J/ipn + n^ 
decay, the xa via its J/V>7 decay, and the r\ c via its K + K~7r° and Ks{t{ + -k~)K~-k + [[| de- 
cays. 

For J ftp and if)(2S) — > decays, we use oppositely charged track pairs where both 
tracks are positively identified as leptons. For the B® — > J /il)Ks{n + ^) mode, the require- 
ment for one of the tracks is relaxed: a track with an ECL energy deposit consistent with 
a minimum ionizing particle is accepted as a muon and a track that satisfies either the 
dE/dx or the ECL shower energy requirements as an electron. For e + e~ pairs, we include 
the four-momentum of every photon detected within 0.05 radians of the original e + or e~ 
direction in the invariant mass calculation. Nevertheless a radiative tail remains and we 
accept pairs in the asymmetric invariant mass interval between — 12.5<r and +3<r of Mj/^ or 
Mm2S), where a = 12 MeV/c 2 is the mass resolution. The fi + radiative tail is smaller; 
we select pairs within —5a and +3cr of Mj/^ or M^s)- Candidate K s — > n + ir~ decays 
are oppositely charged track pairs that have an invariant mass within ±4<r of the K° mass 
(a ~ 4 MeV/c 2 ). For the J/ipKs final state, i^s — > 7r°7r° decays are also used. For 7r°7r° 
candidates, we try all combinations where there are two 77 pairs with an invariant mass 
between 80 and 150 MeV/c 2 , assuming they originate from the center of the run-dependent 
average interaction point (IP). We minimize the sum of the \ 2 values from constrained fits 
of each pair to the 7r° mass with 7 directions determined by varying the decay point along 
the Kg flight path, which is taken as the line from the IP to the energy-weighted center of 
the four showers. We select combinations with a 7r°7r° invariant mass within ~ ±3<r of M^o, 
where a ~ 9.3 MeV/c 2 . For the J/i/jtt° mode, we use a minimum 7 energy of 100 MeV and 
select 77 pairs with an invariant mass within ±3o~ of M n o, where a ~ 4.9 MeV/c 2 . 

We isolate reconstructed B meson decays using the energy difference AE = E^ — E^ m 
and the beam-energy constrained mass Mf, c = \J(E£™ s m ) 2 — (p c j^ s ) 2 , where E£™ s m is the cms 
beam energy, and E c ^ s and p ™ are the cms energy and momentum of the B candidate. 
Figure [I] shows the M^ c distribution for all channels combined (other than J/if)K£) after a 
AE selection that varies from ±25 MeV to ±100 MeV (corresponding to ~ ±3o~), depending 
on the mode. The B meson signal region is defined as 5.270 < Mb c < 5.290 GeV/c 2 ; the 
Mbc resolution is 3.0 MeV/c 2 . Table | lists the numbers of observed events (N ev ) and the 
background (N^gd) determined by extrapolating the event rate in the non-signal AE vs. M^ c 
region into the signal region. 

Candidate B Q d — > J/ipK L decays are selected by requiring the observed K L direction 
to be within 45° from the direction expected for a two-body decay (ignoring the B d cms 
motion). We reduce the background by means of a likelihood quantity that depends on the 
J/ijj cms momentum, the angle between the Kl and its nearest-neighbor charged track, the 
charged track multiplicity, and the kinematics that obtain when the event is reconstructed 
assuming a B + — > J/iI)K* + (Kl7c + ) hypothesis. In addition, we remove events that are 
reconstructed as B° d -> J/tpK s , J/ipK*°(K+-K- , K S 7T°), B + -> J/^+, or J/ijK* + (K + n° 
Ksir + ) decays. Figure |2| shows the pf™ distribution, calculated for a B® — > J/^Kl two- 
body decay hypothesis, for the surviving events. The histograms in the figure are the results 
of a fit to the signal and background distributions, where the shapes are derived from Monte 
Carlo simulations (MC) J7|, and the normalizations are allowed to vary. Among the total of 
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131 entries in the 0.2 < p c ^ s < 0.45 GeV/c signal region, the fit finds 77 J/ipK F events. 

The leptons and charged pions and kaons among the tracks that are not associated with 
fcp are used to identify the flavor of the accompanying B meson. Tracks are selected in 
several categories that distinguish the 6-flavor by the track's charge: high momentum leptons 
from b — * ci~V, lower momentum leptons from c — > s£ + u, charged kaons from b — > c — > s, 
high momentum pions from decays of the type B^ — > D ( -*^(n + , p + , af, etc.), and slow pions 
from D*~ — > D°-k~ . For each track in one of these categories, we use the MC to determine 
the relative probability that it originates from a B® or B& as a function of its charge, cms 
momentum and polar angle, particle-identification probability, and other kinematic and 
event shape quantities. We combine the results from the different track categories (taking 
into account correlations for the case of multiple inputs) to determine a 6-flavor q, where 
q = +1 when f tag is more likely to be a B^ and —1 for a Bd ■ We use the MC to evaluate 
an event-by-event flavor-tagging dilution factor, r, which ranges from r = for no flavor 
discrimination to r = 1 for perfect flavor assignment. We only use r to categorize the event. 
For the CP asymmetry analysis, we use the data to correct for wrong-flavor assignments. 

The probabilities for an incorrect flavor assignment, wi (I = 1,6), are measured directly 
from the data for six r intervals using a sample of exclusively reconstructed, self-tagged 
B® — > D*~i + v, D ( -*^n + , and D*~ p + decays. The 6-flavor of the accompanying B me- 
son is assigned according to the above-described flavor-tagging algorithm, and values of wi 
are determined from the amplitudes of the time-dependent B^-Bf mixing oscillations 0: 
(N OF — N SF )/ (N OF + N SF ) = (1 — 2wi) cos(Am rf At). Here N OF and N SF are the numbers 
of opposite and same flavor events. Table |TI| lists the resulting wi values together with the 
fraction of the events (//) in each r interval. All events in Table | fall in one of the six r 
intervals. The total effective tagging efficiency is I]///(l — 2wi) 2 = 0.270to;o22, wnere t ne 
error includes both statistical and systematic uncertainties, in good agreement with the MC 
result of 0.274. We check for a possible bias in the flavor tagging by measuring the effec- 
tive tagging efficiency for B® and Bd self-tagged samples separately, and for different At 
intervals. We find no statistically significant difference. 

The vertex positions for the fcp and f tag decays are reconstructed using tracks that 
have at least one 3-dimensional coordinate determined from associated rip and z hits in 
the same SVD layer plus one or more additional z hits in other SVD layers. Each vertex 
position is required to be consistent with the IP profile smeared in the rip plane by the 
B meson decay length. (The IP size, determined run-by-run, is typically a x ~ 100 /im, 
o y ~ 5 /im and a z ~ 3 mm.) The fcp vertex is determined using lepton tracks from the 
J/ijj or ip(2S) decays, or prompt tracks from r] c decays. The f tag vertex is determined from 
tracks not assigned to fcp with additional requirements of 5r < 0.5 mm, 5z < 1.8 mm 
and as z < 0.5 mm, where 5r and 5z are the distances of the closest approach to the fcp 
vertex in the rip plane and the z direction, respectively, and o& z is the calculated error of Sz. 
Tracks that form a K$ are removed. The MC indicates that the average zqp resolution is 
75 fim (rms); the z tag resolution is worse (140 /im) because of the lower average momentum 
of the ftag decay products and the smearing caused by secondary tracks from charmed meson 
decays. 

The resolution function R(At) for the proper time interval is parameterized as a sum of 
two Gaussian components: a main component due to the SVD vertex resolution, charmed 
meson lifetimes and the effect of the cms motion of the B mesons, plus a tail component 
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caused by poorly reconstructed tracks. The means (/Vam> Html) an d widths (a main , a tail) of 
the Gaussians are calculated event-by-event from the fcp and ft ag vertex fit error matrices; 
average values are [i main = -0.09 ps, fi ta ii = -0.78 ps and a main = 1.54 ps, a taU = 3.78 ps. 
The negative values of the means are due to secondary tracks from charmed mesons. The 
relative fraction of the main Gaussian is determined to be 0.982 ± 0.013 from a study of 
B® — > D*~l + v events. The reliability of the At determination and R(At) parameterization 
is confirmed by lifetime measurements of the neutral and charged B mesons |J that use the 



same procedures and are in good agreement with the world average values |L0 . 

We determine sin 20! from an unbinned maximum-likelihood fit to the observed At dis- 
tributions. The probability density function (pdf) expected for the signal distribution is 
given by 

-]At|/r s o 

e a. 

V sig (At,q,w h £ f ) = — {1 - £ f q(l - 2w i )sin20isin(Am d At)}, 

2r B o 

where we fix the B® lifetime and mass difference at their world average values ||10|| . The 
pdf used for background events is V k g (At) = f T e~^ At ^ Tbk 9 /2r bkg + (1 — f T )5(At) 1 where 
f T is the fraction of the background component with an effective lifetime r b k g and 5(At) 
is the Dirac delta function. For all fcp modes except J/ipK L we find f T = O.lOtg;^ an d 
Ukg = l-75^o;g2 ps using events in background-dominated regions of AE vs. M bc . The J/tpKL 
background is dominated by B — > J/ipX decays, where some final states are CP eigenstates 
and need special treatment. A MC study shows that the background contribution from the 
£f — —I sources J/ipK s , ip(2S)K s and XciKs is 7.9%, while that from the = +1 ip{2S)K L 
and Xc\Kl modes is 7.0%. Thus, the effects on the CP asymmetry from these states nearly 
cancel. The remaining dominant CP mode, J/iPK*(Kl7t°), which accounts for 19% of the 
total background, is taken to be a 73/27 mixture of = — 1 and +1, respectively, based on 



our measurement of the J/ip polarization in the B® — > J/if)K*°(KsTT°) decay [11]. For the 
J/iPK*(Kltt ) background pdf we use V S i g with effective CP eigenvalue £/ = — 0.46^o!54, 
where the error has been expanded to include all possible values. For the non-CP background 
modes we use Vbkg with f T = 1 and r^g = tb- 

The pdfs are convolved with R(At) to determine the likelihood value for each event as a 
function of sin20i: 

(-'i J {f sigP sigi^At , 

1, wi, 0) + (! _ fsig)V b kg(At')}R(At - At')dAt', 

where f S i g is the probability that the event is signal, calculated as a function of •p B m,s for 
J/ipKi and of AE and for other modes. The most probable sin 2<pi is the value that 
maximizes the likelihood function L = JJiCi, where the product is over all events. We 
performed a blind analysis: the fitting algorithms were developed and finalized using a 
flavor-tagging routine that does not divulge the sign of q. The sign of q was then turned 
on and the application of the fit to all the events listed in Table | produces the result 
sin 20! = 0.58lo!34taio? where the first error is statistical and the second systematic. The 
systematic errors are dominated by the uncertainties in wi (lo!of) an d the J/tl)K L background 
(±0.05). Separate fits to the £/ = -1 and = +1 event samples give 0.82lgj? and O.lOi&io, 
respectively [12]. Figure ^|(a) shows — 21n(L/L 

max ) 

function of sin20i for the = — 1 



and = +1 modes separately and for both modes combined. Figure |3|(b) shows the 
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asymmetry obtained by performing the fit to events in At bins separately, together with a 
curve that represents sin20i sin(AmdAt) for sin20i = 0.58. 

We check for a possible fit bias by applying the same fit to non-CP eigenstate modes: 
B° d -> D*~p+, J/ipK*°(K + ir-), and D*~£ + u, where "sin20i" should be zero, 

and the charged mode B + — > J/ipK + . For all the modes combined we find 0.065 ± 0.075, 
consistent with a null asymmetry. 

We have presented a measurement of the Standard Model CP violation parameter sin 20! 
based on a 10.5 fb" 1 data sample collected at the T(45): 

sin 20x = 0.58^;g(stat)+°;^(syst). 

The probability of observing sin 20i > 0.58 if the true value is zero is 4.9%. Our measurement 



is more precise than the previous measurements [13[ and consistent with SM constraints [JT4 
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TABLE I. The numbers of CP eigenstate events 
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TABLE II. Experimentally determined event fractions (//) and incorrect flavor assignment 
probabilities (wi) for each r interval. 
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FIG. 1. The beam-constrained mass distribution for all decay modes combined (other than 
— > J/ipKi). The shaded area is the estimated background. The dashed lines indicate the 
signal region. 
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FIG. 2. Thep^" s distribution for B® — > J/i/jKl candidates with the results of the fit. The solid 
line is the signal plus background; the shaded area is background only. The dashed lines indicate 
the signal region. 
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FIG. 3. (a) Values of -21n(L/L 

max) vs. sin 2<pi for the £y — —1 and +1 modes separately 
and for both modes combined, (b) The asymmetry obtained from separate fits to each At bin; the 
curve is the result of the global fit (sin20i = 0.58). 
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